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Abstract.  Upper  and  lower  bounds  are  presented  for  the  magnitude  of  the  strain  energy  density 
in  linear  anisotropic  elastic  materials.  One  set  of  bounds  is  given  in  terms  of  the  magnitude  of  the 
stress  held,  another  in  terms  of  the  magnitude  of  the  strain  field.  Explicit  algebraic  formulas  are 
given  for  the  bounds  in  the  case  of  cubic,  transversely  isotropic,  hexagonal  and  tetragonal 
symmetry.  In  the  case  of  orthotropic  symmetry  the  explicit  bounds  depend  upon  the  solution  of 
a  cubic  equation,  and  in  the  case  of  the  monoclinic  and  triclinic  symmetries,  on  the  solution  of 
sixth  order  equations. 


Bounds  on  the  magnitude  ot‘  the  strain  energy  density  in  linear  anisotropic 
elastic  materials  are  needed  in  proofs  of  Saint-Venant’s  principle  for  these 
materials  (see,  for  example,  Toupin  [1],  Horgan  [2],  and  the  review  of  Horgan 
and  Knowles  [3]).  In  this  note  we  extend  the  method  of  bounding  the  strain 
energy  density  employed  by  Horgan  [4]  to  include  specific  anisotropic  elastic 
symmetries.  We  do  this  using  a  result  of  Mehrabadi  and  Cowin  [5]  in  which 
the  coefficients  of  elasticity  are  expressed  as  a  second  rank  tensor  in  a 
six-dimensional  space  rather  than  as  a  fourth  rank  tensor  in  a  three-dimen¬ 
sional  space.  The  results  of  Mehrabadi  and  Cowin  [5]  are  a  development  of 
ideas  due  to  Kelvin  [6];  see  also  Rychlewski  [7],  The  eigenvalues  of  the 
six-dimensional,  second  rank  elasticity  tensor  are  the  numerical  coefficients  in 
the  bounds  obtained. 

The  anisotropic  form  of  Hooke’s  law  is  often  written  in  indicial  notation  as 

Tij  =  Cljkm  Ekm ,  ( 1) 

where  the  C,lkm  are  the  components  of  the  elasticity  tensor.  There  are  three 
important  symmetry  restrictions  on  the  tensor  C,lkm.  These  restrictions,  which 
require  that  components  with  the  subscripts  ijkm ,  jikm,  and  krnij  be  equal, 
follow  from  the  symmetry  of  the  stress  tensor,  the  symmetry  of  the  strain 
tensor,  and  the  requirement  that  no  work  be  produced  by  the  elastic  material 
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in  a  closed  loading  cycle,  respectively.  Written  as  a  linear  transformation  in 
six  dimensions,  Hooke’s  law  (1)  has  the  representation,  T  =  cE,  or 


Tu' 
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in  the  notation  of  Voigt.  The  relationship  of  the  components  of  Ctjkm  to  the 
components  of  the  symmetric  matrix  c  is  easily  seen. 

Introducing  new  notation,  (2)  can  be  rewritten  in  the  form, 

t  =  c£,  (3) 


where  the  shearing  components  of  these  new  six-dimensional  stress  and  strain 
vectors,  denoted  by  T  and  £,  respectively,  are  multiplied  by  Jl,  and  c  is  a 
new  six-by-six  matrix.  Thus  the  matrix  form  of  (3)  is  given  by 
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(4) 


The  symmetric  matrix  c  can  be  shown,  Mehrabadi  and  Cowin  [5],  to  represent 
the  components  of  a  second  rank  tensor  in  a  six-dimensional  space,  whereas 
the  components  of  the  matrix  c  appearing  in  (2)  do  not  form  a  tensor.  The 
inverse  of  the  elasticity  tensor  c  is  the  compliance  tensor  s  where  s  =  c~ ',  thus 


£  =  st. 


(5) 


where  s  is  also  a  second  rank  tensor  in  a  six-dimensional  space.  The  eigenval¬ 
ues  of  the  matrix  c  are  the  six  numbers  A  satisfying  the  equation 


(6) 


(c- Al)Sl  =  0, 

and  their  inverses  are  the  eigenvalues  of  the  matrix  s 
(s  —  ( l/A)l)Sl  =  0, 


(7) 
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where  the  vectors  ft  represent  the  eigenvectors  of  c  or  s.  The  linear  transfor¬ 
mation  (3)  defined  by  c  is  a  six-dimensional  symmetric  transformation, 
assumed  to  be  positive  definite,  which  has,  of  course,  six  positive  eigenvalues. 
These  eigenvalues  will  be  denoted  by  A„  /'=!,.  . ,  6,  and  ordered  by  the 
inequalities  A,  ^  ^  A6  >  0.  It  follows  that 

A6|£|<|f|<A,|£|,  ^-|t|<|£|<  £|f|.  (8) 

where  the  vertical  bars  on  either  side  of  a  vector  indicate  the  norm  of  the 
vector,  e.g. 

|t|  =  /f^t.  (9) 

The  strain  energy  density  is  denoted  by  £  where 

2I  =  t-£  =  (c£)-£  =  T-(st),  (10) 

and  it  can  be  expressed  in  terms  of  the  strain  eigenvector  |£|ft  using  the 
eigenvalues  of  (6)  as 

21  =  |£|2{A,[tf,]2  +  A  2(tf,]2  +  A,[tf3]2  +  A4[tf4]2  +  A5[A5]2  +  AJA,,]2}, 

(ID 

or  in  terms  of  the  stress  eigenvector  jT’jft  using  the  eigenvalues  of  (7)  as 

21 .  |t|.  {1  w  + 1  1*1’  + 1  [*.]■  + 1 1*1’  + 1 1*1’  + 1  [*!’}• 

(12) 

Recalling  that  the  vector  ft  is  a  unit  vector  (ft  •  ft  =  1)  in  six  dimensions,  recall¬ 
ing  also  the  ordering  of  the  eigenvalues  by  the  inequalities  A,  ^  ^  A6  >  0, 

the  results  (11)  and  (12)  yield  the  inequalities 

A6|£|2  ^  2Z  <  A,|£|2,  and  -J-  |f  |2  21  <  ~  |t|2,  (13) 

A|  A6 

respectively.  These  inequalities  represent  the  bounds  of  interest.  For  a  partic¬ 
ular  elastic  symmetry  the  bounds  (13)  are  employed  with  the  values  of  A,  and 
A6  taken  to  be  the  numerically  largest  and  smallest,  respectively,  of  the 
eigenvalues  listed  in  Table  I  for  that  particular  symmetry.  The  second  of  the 
inequalities  (13)  has  been  previously  derived  in  [4]  (p.  232)  and  illustrated 
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Table  l.  The  sets  of  distinct  eigenvalues  A  for  each  of  the  ten  distinct  elastic  symmetries.  The 
multiplicity  of  these  eigenvalues  is  discussed  in  {$]. 

Isotropic  symmetry 

cM  +  2 cl2, 2c*.  (i.e..  32  +  2/i,  2p) 

Cubic  symmetry 

cll  +  C,,  —  C,2,  2c*, 

Transversely  isotropic  symmetry 

i  +  cl2  +  ci3)  +  y8c;,  +  (C|,  +c,2-  c}j)2l. 

^[(c,|  +  C12  +  Cjj)  “  %/8c?j+(C„  +Clj-Cjj)Il 
CI1  ~  CI2>  2c*, 

Hexagonal  (7)  symmetry 
2^ 1 1  +  ^12  +  ^33)  +  n/8cJj  4-  (c ,|  +  c,2  —  Cjj)2]* 

2l(C II  +  C12  +  c3j)  ”  >/^C?3  +  (Cll  +  C12  ~  C33)^l» 

5 Kc  1 1  ““  c\2  +  2^44)  4-  y/\6(c2]4  4-  4*  (C||  — •  cl2  —  2^*44) 

|Kc,,  -  c,2  4-  2c*,)  -  v/|6(C|«  +  C15)  +  (*n  “  ci2  -  2c*,)2]. 

Hexagonal  (6)  symmetry 

2(^1 1  +  CI2  +  c3j)  +  n/8C|j  +  (C,,  +  C,2  —  Cjj)2], 

§KC 1 1  +  CI2  +  Cjj)  —  V/8f  13  +  (CU  +  c12  —  Cjj)2], 

y[(ci  1  -  fu  +  2c*,)  +  >/16ci«  +  (cii  -  C12  -  2c*,)2], 

2l(c<  1  ~  c  1  2  +  2c*,)  -  716c2,  +  (c„  -  c,2  -  2c*,)2], 

Tetragonal  (7)  symmetry 

2l(cll  ■+■  CI2  +  c3j)  +  v/8cfj  4-  (C, |  +C|2  —  C33)2], 
ll(C||  +  C|2  +  C33)  —  \/8Ci3  +  (cll  +  CI2  —  Css)2). 

jKcn  —  C12  +  2 c^*)  +  ^/ibc2^  +(c„  - c,2  —  2CJJ)2], 

2((cu  —  ci2  +  2cw  —  ^16c26  +  (c  1 1  —  c,2  —  2c,*)'),  2c*, 

Tetragonal  (6)  symmetry 

i((c„  +  C|2  +  C3j)  +  V8c?j  +(C|,  +  C,2  —  Cjj)‘|, 

2l(C||  +  C|2  +  Cjj)  —  >/8c2j  +  (c,,  +  C|2  —  Cjj)'J, 

C||  -c,2.  2 cM,  2c*, 
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Table  I  ( Continued ) 


Orthotropic  symmetry 


2c*.  2c5J,  2c66,  and  the  eigenvalues  of  the  main* 


Monoclinic  symmetry 


Eigenvalues  of  the  matrix 


Triclinic  symmetry 


Eigenvalues  of  the  matrix 
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there  for  the  case  of  isotropic  symmetry.  The  first  of  (13)  has  been  given  in 
Gurtin  [8]  (p.  85)  and,  in  the  case  of  isotropic  symmetry,  by  Villaggio  [9]  (p. 
46).  Toupin  [1]  and  Gurtin  [8]  employ  the  terminology  maximum  elastic 
modulus  for  A,  and  minimum  elastic  modulus  for  A6;  see  also  [3]  (p.  240). 
Kelvin  [6]  called  the  A„  i  =  1 , . . . ,  6,  the  six  principal  elasticities  of  the 
material;  Pipkin  [10]  uses  the  term  “principal  compliance”  for  the  inverse  of 
the  same  quantities;  and  Rychlewski  [7]  suggests,  with  persuasive  historical 
justification  (but  contrary  to  the  contemporary  trend  to  avoid  eponyms),  that 
they  be  called  the  Kelvin  moduli.  Applications  of  the  inequality  represented  by 
the  first  of  ( 13),  in  the  isotropic  case,  to  obtain  bounds  for  total  energies  and 
related  quantities  are  described  by  Villaggio  [9]  (see,  e.g.,  p.  400.418). 

In  closing,  we  remark  that  the  positivity  of  i  =  I,., . . ,  6,  is  equivalent  to 
the  positive  definiteness  of  £,  expressed  as  a  quadratic  form  in  either  "f  or  £. 
The  explicitness  of  the  eigenvalues  tabulated  in  Table  1  here  thus  enables  one 
to  write  down  necessary  and  sufficient  conditions  for  positive  definiteness  of  £ 
directly  in  terms  of  the  symmetric  Voigt  matrix  c.  The  resulting  conditions  arc 
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simpler  than  the  expressions  involving  principal  minors  of  c  that  are  usually 
employed  in  the  literature. 
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